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Abstract 

The quantization of the scalar and electromagnetic fields in the presence of a parabolic 
mirror is further developed in the context of a geometric optics approximation. We extend 
results in a previous paper to more general geometries, and also correct an error in one 
section of that paper. We calculate the mean squared scalar and electric fields near the 
focal line of a parabolic cylindrical mirror. These quantities are found to grow as inverse 
powers of the distance from the focus. We give a combination of analytic and numerical 
results for the mean squared fields. In particular, we find that the mean squared electric 
field can be either negative or positive, depending upon the choice of parameters. The case 
of a negative mean squared electric field corresponds to a repulsive Van der Waals force on 
an atom near the focus, and to a region of negative energy density. Similarly, a positive 
value corresponds to an attractive force and a possibility of atom trapping in the vicinity 
of the focus. 
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I. INTRODUCTION 



In a previous paper henceforth I, we developed a geometric optics approach 
to the quantization of scalar and electromagnetic fields near the focus of a parabolic 
mirror. We found that there can be enhanced fluctuations near the focus in the 
sense that mean squared field quantities scale as an inverse power of the distance 
from the focus, rather than an inverse power of the distance from the mirror. These 
enhanced fluctuations were found to arise from an interference term between different 
reflected rays. In the present paper we extend the previous treatment. In I, only 
points on the symmetry axis were considered. Here we are able to treat points in an 
arbitrary direction from the focal line of a parabolic cylinder. We give more detailed 
numerical results which provide a fuller picture of the phenomenon of focusing of 
vacuum fluctuations. We also correct some erroneous results in Sect. V of I. 

The outline of the present paper is as follows: In Sect. ||, we review and extend 
some of the formalism used to calculate the mean squared scalar and electric fields, 
{ip^) and (E^), near the focus. In Sect. |T| we derive some geometric expressions 
which are needed to study fluctuations at points off of the symmetry axis. In Sect. IV, 
the evaluation of integrals with singular integrands is revisited. Two different, but 
equivalent, approaches are discussed. A particular case where the integrals can be 
performed analytically is examined in Sect. 0. More generally, it is necessary to 



calculate (yj ) and (E ) numerically. A procedure for doing so is outlined in Sect. ^ 



Some detailed numerical results are also presented there. The limits of validity of our 
model and results will be examined in Sect. |V11| . This discussion will draw on some 
results on diffraction obtained in the Appendix. The experimental testability of our 
conclusions will be discussed in Sect. |VIII| . Finally, the results of the paper will be 
summarized in Sect. 

Units in which % = c = 1 will be used throughout this paper. Electromagnetic 
quantities will be in Lorentz-Heaviside units. 



II. BASIC FORMALISM 

Here we will briefly review the geometric optics approach developed in I. The basic 
assumption is that we may use a ray tracing method to determine the functional form 
of the high frequency modes, which will in turn give the dominant contribution to 
the expectation values of squared field operators. We start with a basis of plane wave 
modes. The incident wave, for a scalar field, may be taken to be 

with box normalization in a volume V . In the presence of a boundary, this is replaced 
by the sum of incident and reflected waves, 

i^k = ^ + E/k^ (2) 
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where the fj^^ are the various reflected waves. One could also in principle adopt a 
wavepacket basis, in which Fk is replaced by a localized wavepacket. Because the time 
evolution preserves the Klein-Gordon norm, if the various modes are orthonormal in 
the past, they will remain so after reflection from the mirror. Thus we can view 
Eq. as the limit of a set of orthonormal wavepacket modes in which the modes 
become sharply peaked in frequency and hence delocalized. 

It was shown in I that the renormalized expectation value of the squared scalar 
field is given by a sum of interference terms: 



i i^j 



(3) 



This renormalized expectation value is defined as a difference in the mean value of 
(/)^ with and without the mirror, and hence will vanish at large distances from the 
mirror. The various interference terms in the above expression yield contributions to 
{ip'^) which are proportional to the inverse square of the appropriate path difference. 
Thus in the vicinity of the focus, the interference terms between different reflected 
rays will dominate over that between the incident and a reflected ray. In the present 
paper, we will consider cases with no more than two reflected rays, and write 

{^') ^2ReY. fi'^fK. (4) 

k 

In the case of a parabolic cylinder, this may be expressed as 

Here Ai is the path difference for the two reflected rays, and the integration is over 
the reflection angle of one of the rays. The range is chosen so that each pair of 
reflected rays is counted once. The corresponding expression for the mean squared 
electric field near the focus of a parabolic cylinder is found in I to be 



III. OPTICS OF PARABOLIC MIRRORS 

In this section, we wish to generalize some of the results of I concerning the incident 
and reflected rays in the presence of a parabolic mirror. Consider the geometry 
illustrated in Fig. |l|. An incident ray at an angle of 6 to the symmetry axis is 
reflected from the point {xi,yi), and then reaches the point P at an angle of 6'. We 
first need to find the relation between 6 and 6'. Note that the reflected ray crosses 
the symmetry axis at a distance c from the focus. It was shown in I that 

b{l-cose')' ^' 
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FIG. 1: An incident ray at an angle 6 reflects off of a parabolic mirror at an angle 9' and 



arrives at the point of interest P. This point is a distance a from the focus F in a direction 
at an angle 7 with respect to the symmetry axis. Throughout this paper, we assume that 
a <C 6, where is the distance from the focus to the mirror. 

However, from the law of sines, we have that 

c sin 6*' = a sm{9' — 7) . (8) 

Hence we now have that 

o = lm, (9) 

where 

There will be multiply reflected rays whenever different values of 6' are associated 
with the same value of /. The function f{6') is plotted in Fig. ^ for various values 
of 7. We can see from these plots that in general there can be up to four reflected 
angles 6' for a given incident angle 6. However, if the mirror size ^0 is restricted to 
be less than 27r/3, then there will never be more than two values of 6' for a given 
6. Throughout this paper, we will assume 6q < 27r/3, and hence have at most two 
reflected rays for a given incident ray. The two reflected rays will occur at ^' = a 
and 9' = /?, where 

/(«) = /(/?). (11) 
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FIG. 2: The function f{9') is plotted for various values of 7. This function relates the 
angle of the incident ray, 9, to the angle of the reflected ray, 9', through the relation 
9 = {a/b)f{9'). 



Our next task is to calculate the difference in path length, Ai, for these two 
reflected rays. Again, this is a generalization of a calculation given in I. Consider the 
situation illustrated in Fig. ^, where a reflected ray with angle 6' reaches the point P 
after reflecting from the point (xj, t/i) on the mirror. Let Si be the distance traveled 
after reflection, and S2 be the distance traveled between when the ray crosses the line 
X = Xo and when it reaches the mirror. Note that 

S2 = {xi - Xo) sec6 ^ Xi- xo + . (12) 

Even if Xi < 0, we can choose xq to be such that S2 > 0. Next note that 

Si = {i/i — asin.'y) csc9' . (13) 
The reflection point {xi,yi) is the intersection of the line 

y = tan6''(a; — acos7) + asin7 (14) 
with the parabola, given by 

X = — - — . (15) 
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FIG. 3: The distance travelled by a ray which reaches the point P is illustrated. Here 
the line x = xq is an arbitrarily chosen vertical line somewhere to the left of the point of 
reflection. The incident ray travels a distance S2 from this line to the point of reflection. 
After reflection it travels a distance si before arriving at the point P. This point is a 
distance a from the focus F in a direction given by the angle 7. 



These relations lead to 



Vi — b cot 9' 



sec e\j 1 - 2 ( ^ ) sin^ ^'(cos 7 - sin 7 cot 9') - 1 



(16) 



We may now expand this expression to first order in a/6 and combine it with our 
previous expressions to show that, to first order, 



I — Si -\- S2 — h — Xq — a(cos 7 cos 9' — sin 7 sin 9') . 



(17) 



Thus the magnitude of the path length difference for a and /3, two different values of 
9', is 

/S.l — a\ cos7(cosq; — cos/3) + sin7(sinQ; — sin /3)| . (18) 
In the hmit that 7 = 0, we obtain the result for used in I. 
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IV. EVALUATION OF SINGULAR INTEGRALS 



We can express Eq. (1^) as 



and Eq. @ as 



1 



da 
da 



(19) 



(20) 



Here 



h = h{a, P) = I cos 7(cos a — cos P) + sin 7(sin a — sin /3) | , (21) 



and the integrations run over the full range for which there are two reflected rays, 
thus counting each pair twice. 

The integrands of these integrals contain singularities within the range of integra- 
tion. These singularities occur at a critical angle, 6c, at which h vanishes linearly. 
Thus the (y?^) integral contains a (a — 6'c)~^ singularity, and the (E^) contain a 
(a — 6'c)~'' singularity. These singularities are presumably artifacts of assuming per- 



fect mirrors with sharp boundaries as will be discussed in Sect. [VII[ The singularity 
occurs when both a and (3 are approaching 6c from opposite directions. The critical 
angles 6c are just the extrema of the function f{6'). 

Integrals with such singular integrands can be defined by a generalization of the 
principal value prescription. This generalization involves an integration by parts to 
recast the original integral as one containing a less singular integrand, plus surface 
terms. For example, we may use 



2 dx' 



Inx 



to write 



['dx^ = t dxf\x) lna;2 + 

J a X^ 2 J a 



+ o/ (^) l'^^ 

X 2 



Similarly, 



leads to 



x^ 



r-r In , 

12 dx* 



(22) 



(23) 



(24) 



x'' 



6X'^ ox^ bx 12 



12 Ja dx^ 



(25) 

If 6 > and a < 0, there were singularities in the original integrals which are replaced 
by integrable, logarithmic singularities. In all cases, the surface terms are evaluated 
away from the singularity and are hence finite. Note that we could have added 
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additional polynomial terms on the righthand sides of both Eqs. (^) and (|^). For 
example, we could replace hix^ by Inx'^ + cix + co in Eq. (p^ . However, the arbitrary 
constants cq and ci will cancel out between the two terms on the righthand side of 
Eq. p3| ) and hence can be ignored. A similar cancellation occurs if we add additional 
terms into Eq. (P^). 

If / and a sufficient number of its derivatives vanish at the endpoints, then the 
surface terms vanish. In I, it was incorrectly argued that the surface terms which 
arise in the present problem can be ignored. It is true that if the reflectivity of the 
mirror falls smoothly to zero, the surface terms vanish. However, in this case it would 
be necessary to integrate explicitly the function which is falling smoothly to zero. If 
the reflectivity falls very rapidly at the edge of the mirror, its derivatives will be large, 
and effectively reproduce the surface terms. Thus, the explicit results given in Sect. 5 
of I are incorrect. In Sects. and |V| of the present paper, we will give new results 



which replace and generalize the older results. 

There is an alternative way to implement the above integration by parts prescrip- 
tion. This is simply to evaluate the integral as an indefinite integral, and evaluate the 
result at the endpoints, ignoring the singularity within the integration range. Even if 
it is not possible to find a closed form expression for the indefinite integral, one can 
expand the integrand in a Laurent series about x = and integrate term by term, 
using the relation 

''^d^ = -^(-^-i:^) ■ (26) 



n + 1 Va"+i 

One can verify this relation for the case a < < 6 using the above integration by 
parts method. 



The integrands in Eqs. (19) and (^) arise from the integrals 

duu cosux (27) 



oo 



1 

Jo 
and 

1 1 

— = - / duj'^ uj coscjx, (28) 
x^ 6 JO 

respectively. The integrals are understood to be performed with the aid of a con- 
vergence factor, such as e~'*"^, with the limit A ^ taken after integration. The 
singularities at x = are due to the contributions of arbitrarily large values of u. If 
there were to be a physical cutoff at high frequencies, then these singulariries would 
disappear. However, the integrals containing the singular integrands are independent 
of the cutoff in the limit that it occurs at sufficiently high frequencies. As an example, 
consider the function 

g2{x,X)= duuje'^'^ cos ux = —-— — r-^- . (29) 
Jo (A^ + x^j^ 

As A — s> 0, g2{x, A) — but for A 7^ 0, g2{x, A) is finite for all x. Furthermore, 
if we integrate g2{x, A) on x through x = 0, and then take the A — > limit, the result 
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will be finite and the same as that obtained by the above formal procedures: 



XQ 2xo 2 

92{x,X)dx= I ^2 ^ — • (30) 

-xo A -\- Xq Xq 



Similarly, if we define 



g,{x,\)=J^ duu'e-^'^ cosiJx = ^ (PTW ' ^ ^ 

then gi{x, A) as A ^ 0. However, 

g,ix,X)dx = -—-3 .--3. (32) 

Thus we can understand the physical meaning of the formal integration procedures 
discussed above as follows: they provide a shortcut method to obtain the results one 
would obtain by inserting a frequency cutoff in the u integrals, and then removing 
this cutoff after integration over x. 

One peculiar feature of our results will be that both (ip'^) and (E^) diverge for 
particular values of 7. This arises when 6c, the extremum of f{0'), approaches the 
edge of the mirror, 6' = ±6q. The mathematical reason for the divergence is that 
one limit of integration is approaching a point at which h = and the integrand is 
singular. This corresponds to letting either a — with b fixed, or 6 — with a 
fixed, in Eq. (|26|). The physical origin of this singular behavior is that we have made 
two unrealistic assumptions. The mirror is assumed to be perfectly reflecting at all 
frequencies and to have sharp edges at 9' = ±9q. This will be discussed further in 
Sect. 



V. EXACT RESULTS FOR 7 = f 

In general, the integrals for {ip"^) and (E^), Eqs. (|^) and (H), respectively, can only 
be evaluated numerically. This is in part because the second reflection angle (3 is 
implicitly given as a function of the first angle, a. There is one case in which this 
relation may be written down in closed form. This is when 7 = ^, so the function 
f{9') is symmetrical about the origin. In this case, we have 



and we can write 



p = -a, (33) 



2\ 1 da cot^o .O.N 

= 3 2 / — ^ = 103 2 • 34 



Similarly, 



/E^\ = - f^" d^ = cos go(3-2cos^ 9o) 
^ ' 207r3 a4 J-eo sin^ a 307r3 sin^ 9o ' ^ ' 
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FIG. 4: The exact solutions for (99^) and (E^) for the case that 7 = f are plotted as 
functions of the mirror size. Note that a^{ip^) and a^(E2) are dimensionless in our units. 

These expressions are plotted in Fig. |. Note that both {ip^) and (E^) can have 
either sign: 

((^2) > and (E^) < for < < | , 
{ip") < and (E^) > for ^ < ^0 < ^ • (36) 
Both quantities seem to diverge in the limit that ^0 ^ 



127r3 a2 Oo 



This divergence is a special case of the singularities arising when the edge of the 
mirror approaches an extremum of f{0'). 

We can go further and find the derivatives of both (ip^) and (E^) with respect to 
7, at 7 = |. First, we need to find /3(a) to first order in 7 — |. We may expand the 
relation /(/?) = /(a) to this order and show that 

2sin2a(7-^) , , 

(5 = -a- -j^ ^ + ■ ■ ■ . 38 

(cos a — 1) (2 cos a + 1) 
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We can next expand 1/h? and l/h'^ to first order in 7 — ^ . The first order corrections 
are odd functions of a whose exphcit forms we will not need. 

The crucial effect of varying 7 slightly away from 7r/2 is to change the range of 
integration. First consider the case where 7 < |. The integration range now becomes 



where 



- 9q < a < 9q - 5 , 

2sin^ go(f -7) 
;i-cos^o)(2cos^o + l' 



(39) 
(40) 



To first order in 7 — ^, we can write 



r^o-s da 



1 



do 



da 



do 



+ 



da 



00 



do 



da 



00-5 



(41) 



where the subscripts on l//i^ denote the order in an expansion in powers of 7 — ^. 
The first term on the right hand side of Eq. (^) is just the zeroth order part given 
by Eq. (0). The next term vanishes because {l/h?)i is an odd function. The final 
term may be approximated using 



00 



da 



00-S 



1 



This leads to 



127r3 a2 



S 



-7 



tan 6*0 (1 — cos6'o)(2 cos^^o + 1' 



0=6*0 



+ • 



7 < 



TT 



(42) 



(43) 



We may now repeat this procedure for 7 > f . In this case, the range of integration 



IS 



- ^0 - ^ < tt < ^0 , (44) 
and 5 < 0, so the range of integration has again decreased. This causes the first order 
change in (yj^) to have the same magnitude but the opposite sign from the previous 
case. Thus in all cases, we can write 



1 



127r3( 



17- tl 



tan 6*0 (1 — cos 6^0) (2 cos 6*0 + r 



+ ■ 



(45) 



A similar analysis may be applied to (E^) with the result 



57r3 



cos 6^0(3 — 2 cos^ 9o 
24 sin^ Oa 



17-^1 



8sin^ ^0(1 - cos ^o) (2 cos ^0 + 1^ 



+ 



(46) 

Thus both (y)^) and (E^) have cusps at 7 = |. The nonanalytic behavior is due 
to the fact that the range of integration decreases whenever 7 moves away from | 
in either direction. As will be discussed in Sect. |Vli| , the cusp is presumably an 
artifact of a sharp edge approximation, and should be smoothed out in a more exact 
treatment. 
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VI. NUMERICAL PROCEDURES AND RESULTS 



Apart from the special cases discussed in the previous section, it is necessary to 
evaluate {tp^) and (E^) numerically. The first step is to find the second reflection 
angle /? as a function of the first reflection angle a. This involves a straightforward 
numerical solution of the equation 

= /(«) • (47) 
Because we assume a restriction on the angle size of the mirror, 

^0 < y , (48) 

there will always be either one or no roots for (3. Within the geometric optics ap- 
proximation that we use, integrands are assumed to vanish in regions where there are 
no roots. 

As was discussed in Sect. there are at least two methods that may be used for 
explicit evaluation of the integrals which appear in the geometric optics expressions 
for {if^) and (E^). The first is an integration by parts, which replaces the singularity 
in Eqs. (^) and (|20| ) by a logarithmic singularity. Here we will outline how this may 
be done explicitly. Consider first Eq. (|19|), which may be expressed as 

^ " ~6^I [h{a,P{a))]^ ' ^^^^ 

Suppose that we are interested in integrating over the range ai < a < a2 and that 
this range contains one point a = ao at which h = 0. We can write 

/ = io 2 2 / da — — \og{a - ao) ■ (50) 



127r2 Jai [h{a,y{a))]'^ da^ 
The quantity (a — a^Y l[h{a, y{a))f is finite at a = ao- This expression may now be 



integrated by parts, as in Eq. (^3]) . The derivatives of h which arise are computed 
using Eq. (^) and the relation 

da /'(/?)• ^ ^ 

An analogous procedure can be applied to (E^). The detailed expressions for ((^^) and 
(E^) which result from this procedure are rather complicated, and will not written 
down explicitly. 

The second method which may be employed is a variant of the direct integration 
illustrated in Eq. (pG]). The actual integrands in Eqs. (|19D and ( pO]) are too complex 
to integrate in closed form. However, we can use numerical integration of l/h? and 
in regions away from zeros of h and direct integration of a series expansion in 



12 



the neighborhood of a zero. The first step in the generation of the series expansion 
is to expand P{a) around a = ao, 



f3 = 2aQ — a + a2{a — ao) + a^{a — ao) +■■■• (52) 

The coefficients are found by inserting this expansion into Eq. ( ^Tj) and then expand- 
ing both sides of the resulting expression in powers of a — a^. A few of the leading 
coefficients are 



02 
Os 



f" 



3/" 

irr 



where the derivatives of / are evaluated at a = Oq- The expansion for y{a) is then 
used to generate analogous expansions for l/h"^ and which are in turn integrated 
over the interval ~ < « < «o + using Eq. (^61). The result is combined with 
the direct numerical integration outside of this interval. Here is an arbitrary small 
positive number. One of the tests of the numerical procedure is the independence of 
the results upon the choice of ^o- 

We have developed numerical routines based upon both of the above procedures. 
The plots which are given below were created using a routine based upon the second 
method, with an expansion of y to sixth order in a — ao- The particular values of 
^0 were in the range 0.1 < < 0.3. In this range, the routine is relatively stable 
and insensitive to .^o- Smaller values of .^o can lead to instability, as the contributions 
from |x — XqI > 'Co and from |a; — Xo| < are both large in magnitude and tending 
to cancel each other. 

Numerical results for and (E^) for various values of are shown in Figs. |-12 



as functions of 7. In all cases, there are values of 7 at which or (E^) are singular 
in the model of a perfectly reflecting mirror with sharp edges. These points occur 
when an extremum of the function f{6) sits at an edge of the mirror, 6' = ±6o- The 
values and slopes of (ip"^) and (E^) at 7 = | are in agreement with Eqs. (^5]) and 
(^). We can see from the graphs that for smaller mirrors, 6q < ^, we have {(f^) > 
and (E^) < everywhere. For larger mirrors, 60 > ^, we have regions where {(p"^) > 
and other regions where it is negative, and similarly for (E^). The signs of {(f"^) 
and of (E^) always seem to be opposite. The regions near 7 = f for larger mirrors 
where (E^) > are of special interest. These are regions where an atom will feel an 
attractive force toward the focus, and thus has the possibihty of being trapped. 
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FIG. 7: {^'^) for = 1.8. 
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FIG. 8: {ip^) for 9q = 1.8, showing in detail the region where {ip^) < 0. 
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FIG. 11: (E2) for 6*0 = 1.8. 
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FIG. 12: (E^) for = 1-8, showing in detail the region where (E^) > 0. 
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VII. LIMITS OF VALIDITY OF THE RESULTS 



In this section, we will discuss the likely ranges of validity of the model we have 
used to calculate {ip^) and (E^). In particular, we have assumed a geometric optics 
approximation, and a mirror which is both perfectly reflecting and has sharp edges. 
Each of these assumptions will be examined critically. 



A. Geometric Optics 

The use of geometric optics amounts to ignoring diffraction effects, so we can 
gauge the accuracy of geometric optics by estimating the size of these effects. In the 
geometric optics approximation, the incident and reflected waves in Eq. (H) all have 
the same magnitude. Suppose that we now introduce a correction due to diffraction 
and write 

Fk = /k + E/f^ + ^^k. (54) 

i 

In the case of a plane strip which has width 2?/o in one direction and is infinite in 
the other direction, the magnitude of the diffraction correction is estimated in the 



Appendix with the result |AFk/Fk| ^ y^/vo^ where A is the wavelength of the 
incident wave. In our case, we have a parabolic cylinder characterized by the length 
scale b. We could imagine representing the parabolic cylinder by a set of strips which 
are infinite in the z-direction. The order of magnitude of the diffraction correction 
should be the same as for a single strip of width b. Thus we estimate that in the 
present case 




(55) 

The wavelengths which give the dominant contribution to (y^^) and (E^) near the 
focus are those of order a. Thus the interference term between AF^ and any of the 
other terms in Eq. (^) should yield contributions to (y^^) and (E^) which are smaller 

than the dominant contribution by a factor of the order of ^a/b. Thus we estimate 
the diffraction contribution to (ip"^) to be of order 

A(^')~^, (56) 



a2 62 



and that to (E^) to be of order 



A{E')^^. (57) 

aa 02 



Note that the geometric optics results near the focus depend only upon a and 60, 
the angular size of the mirror, but not upon b, the linear dimension of the mirror. 
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However, the diffraction correction decreases with increasing b and hence can be made 
smaller for a larger, more distant mirror. 

In classical optics, the effects of diffraction are normally of the order of the wave- 
length divided by the size of the object. This arises when one is looking at the power 
in the diffracted wave, which is given by the square of Eq. In our case, there 

is a possible contribution from an interference term between the diffracted wave and 
the geometric optics contributions. If a more detailed calculation were to find that 
this term vanished, then our estimates for A(y9^) and A(E^) would become smaller 

than Eqs. ( ^6|) and (|57|) , respectively, by an additional factor of Ja/b. 



B. Edge Effects and Finite Reflectivity 

Our model of the mirror is one in which it is not only perfectly reflecting, but also 
has a sharp edge at which the reflectivity falls discontinuously to zero. Both of these 
assumption are over simplifications to which the singular behavior of {ip"^) and (E^) 
may be attributed. Consider the sharp edge assumption. A more realistic model 
might have the reflectivity falling smoothly to zero over an angular interval of width 
A6 at the edges of the mirror. This would remove the singularities found above at 
specific values of 7. Recall that these singularities arise when an edge of the mirror, 
at ±^0) sits at an extremum of f{0'). In the case of a smoothed edge, both (ip^) and 
(E^) will be bounded for fixed a 

l(.^>l < ^ (58) 

and 

Smoothing of the edges of the mirror is also expected to remove the cusps at 7 = ^. 

The smoothed edges do not, however, remove the singularities as a — > 0. This sin- 
gularity is presumably due to the assumption of perfect reflectivity at all wavelengths. 
A more realistic model would have the reflectivity go to zero at short wavelengths. 
Suppose that the mirror becomes transparent for wavelengths less than some mini- 
mum value, \m- Then we expect to find the bounds 

l(^')l<^ (60) 



m 



and 



l(E^)l<^- (61) 



One reason for reduced reflectivity at short wavelengths is dispersion. The mirror can 
be regarded as close to perfectly reflecting only for wavelengths longer than about 
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the plasma wavelength Ap of the metal in the mirror. Thus our results are only valid 
when a > Xp. For aluminum, for example, Ap ~ 84nm. However, it is unlikely that 
dispersion alone is capable of removing all short distance singularities. The reason 
for this is that dielectric functions approach unity as (A/Ap)^ as A — 0. This is not 
fast enough to regulate integrals which diverge quartically as the short wavelength 
cutoff. In the case of a plane interface between vacuum and a dispersive material, 
quantities such as (E^) still diverge as the interface is approached, although less 
rapidly than in the case of a perfect mirror 0. Thus it seem that other effects, 
such as surface roughness, the breakdown of the continuum description at the atomic 
level, or quantum uncertainty in the location of the interface [|] are needed to produce 
finite values of (E^). If any effect does produce a sufficiently sharp cutoff at short 



wavelengths below about Am, then we would have the bounds in Eqs. ( |60D and (|6l|) 
even without removing the assumption of a sharp edge to the mirror. 



VIII. OBSERVABLE CONSEQUENCES 

In I, several possible experimental tests of the enhanced vacuum fluctuation near 
the focus were discussed. Here these tests will be reviewed and further discussed. The 
most direct test would seem to be to measure the force on an atom or other polarizable 
particle. In a regime where the atom can be described by a static polarizability a, 
the force is F = —W, where 

V = -^a{E'). (62) 
In the vicinity of the focus, we have found that 

(E^) = ^ , (63) 

where A is a dimensionless constant. For an atom in its ground state, Eq. (|62|) should 
be a good approximation when a is larger than the wavelength associated with the 
transition to the first excited state. 

One might try to measure the defection of atoms moving parallel to the focal line 
of a parabolic cylinder. The analogous experiment for a fiat plate was performed by 
Sukenik et al |Q and confirmed Casimir and Polder's theoretical prediction 0. In 
the present case, the expected angular deflection is 

Here m^va = 3.8 x 10^^'^gm and ajva = 3.0 x lO^^^cm^ denote the mass and polar- 
izability of the sodium atom, respectively. (Note that polarizability in the Lorentz- 
Heaviside which we use is 47r times that in Gaussian units.) If t is of order 10~^s 
(the time needed for an atom with a kinetic energy of order 300K to travel a few 
centimeters), and z is of order l/xm, the fractional deflection is significant. 
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An alternative to measuring the deflection of the beam might be to measure the 
relative phase shift in an atom interferometer, which are sensitive to phase shifts of 
the order of 10~^ radians 0]. If one path of the interferometer were to be parallel to 
the focal line at a mean distance of a and the other path far away, the phase difference 
will be 

Both deflection and phase shift measurements would require the atoms to be rather 
well coUimated and localized near the focal line. 

Perhaps the most dramatic confirmation of enhanced vacuum fluctuations would 
be the trapping of atoms near the focus. This would require a mirror with > f so 
that there is a region with (E^) > 0. It would also require the atoms to be cooled 
below a temperature of about 

(i^) (^^) ■ (-) 

Thus at temperatures of the order of 10~^i^, atoms could become trapped in a region 
of the order of l/im from the focus. This type of trapping would be quite differ- 
ent from that currently employed [0 in that it would require no applied classical 
electromagnetic fields. 

The peculiar property that the force on an atom can be attractive from certain 
directions and repulsive from other directions requires some discussion. If one ap- 
proaches from an attractive direction, the potential energy is becoming increasingly 
negative, V ~ — whereas from a repulsive direction it is becoming large and 
positive, V ~ +l/a^. In both cases, ^ ^ for large a. At some minimum value of 
a, the inverse fourth power behavior has to be modified. Thus the exact potential 
energy V{a, 7) must be a continuous function with a stable minimum, as illustrated 
in Fig. [T^. Note that the requirement that — > for large a in all directions rules 
out the possibility that there is only a saddle point, and no minimum. Although 
we can infer the existence of the minimum from our results, the detailed form of 
V{a, 7) near the minimum could only be found by an analysis which goes beyond the 
approximations made in the present paper. 



IX. DISCUSSION AND CONCLUSIONS 

In this paper we have further developed a geometric optics approach to the study 
of vacuum fluctuations near the focus of parabolic mirrors. The main result of this 
approach is that the mean squared scalar field (y?^) and mean squared electric field 
(E^) grow as inverse powers of a, the distance from the focus, for small a. The key 
justification of geometric optics is its self-consistency. When {lp^) and (E^) are large, 
the dominant contributions must come from short wavelengths for which geometric 



optics is a good approximation. This was discussed more quantitatively in Sect. [VII 
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FIG. 13: A possible form of the potential energy function V{a, 7) is illustrated. The solid 
hnes show the region where V can be rehably computed within our approximations. The 
soHd line on the right is the form of y as a function of a for a value of 7 for which the 
force is attractive. The soHd hne on the left is for a value (not necessarily in the opposite 
direction) for which it is repulsive. The dashed line is a possible form which interpolates 
between these two regions, but which would have to be computed by a more detailed theory. 

We have given some explicit analytic and numerical results for {(f"^) and (E^) 
near the focal line of a parabolic cylindrical mirror. In this paper, we restricted our 
attention to the case that the angular size of the mirror, 29q, is less than This 
insures that there are never more than two reflected rays for a given incident ray 
and simplifies the analysis We find that for smaller mirrors, 6q < ^, {ip'^) > and 
(E^) < everywhere. In this case, an atom will feel a repulsive force away from the 
focus. For larger mirrors, | < 6*0 < these mean squared quantities can have either 
sign, depending upon the direction from the focus. In directions nearly perpendicular 
to the symmetry axis, 7 ~ |, we find (y?^) < and (E^) > 0. In this case, the force 
on an atom is attractive toward the focus and trapping becomes a possibility. 

In the geometric optics approximation, the mean squared electric and magnetic 
fields are equal, so the local energy density is equal to (E^). Thus, when (E^) < 0, 
the local energy density is negative, and one has columns of negative energy density 
running parallel to the focal line of the parabolic cylinder. 

The feasibility of experiments to observe the effects on atoms near the focus was 
discussed in Sect. [VII]| . Although the effects are small, it seems plausible that they 
could be observed. 
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APPENDIX 

In this Appendix, we will develop a method which goes beyond the geometric optics 
approximation, and use it to estimate the size of the corrections due to diffraction 
effects. Here we will discuss only a scalar field which satisfies Dirichlet boundary 
conditions on the mirror. The basic method is to write down and then approximately 
solve an integral equation for the scattered wave. This method was first used by 
Kirchhoff to discuss diffraction in classical optics. Let 0(x) be a solution of the 
Helmholtz equation 

V20 + fc20 = O, (A.l) 
and let G(x, x') be a Green's function for this equation 

V'G + A;' G = 5(x - x') . (A.2) 

If we multiply Eq. ( [A.l| ) by G, Eq. ( |A.2| ) by 0, take the difference, and integrate over 
an arbitrary spatial volume V, the result is 

0(x') = f da - [0(x) VG(x, x') - G(x, x') V0(x)] . (A.3) 

s 

Here S is the boundary of V, and da is outward directed. 

We are interested in calculating the scattered wave from a open surface 5*1, rather 
than the solution on the interior of a closed surface S. However, we can let the 



closed surface 5* consist of three segments, as illustrated in Fig. 0. The first is the 
surface of interest 5*1, the second is a segment 5*2 which is hidden from the view of 
an observer in the region where we wish to find the scattered wave, and the third 
closes the surface at a large distance. We now ignore the contributions of 5*2 and S3, 
and assume that = on 6*1, so that we have 



/ rfa ■ V0(x) G(x, x') . (A.4) 

J Si 



iSi 

This is still an integral equation which relates values of (p on the surface to those 
off of the surface. We can solve it in an approximation in which an incident wave 
scatters only once from the surface. Let 0j„ be the incident wave, and assume that 
on Si 

n ■ V0 = 2n ■ V(j)in ■ (A.5) 
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FIG. 14: A closed surface S consists of three parts, the open surface of interest 5i, a 
segment 5*2 hidden from the view of an observer at x', and a portion 5*3 at a very large 
distance. An incident wave (pin reflects from Si and produces the scattered wave (ps- Here 
X is an arbitrary point in ^i. 

We next choose the empty space Green's function 

ifc|x-x'| 

G'(x,x') = --^ -, (A.6) 

47r|x — X I 

and interpret Eq. (|A.4|) as giving the scattered wave, 0^, 

-2 / rfa- V0,„(x)G(x,x'). (A.7) 



fx') 



The physical interpretation of this expression is that each point on 5*1 radiates a scat- 
tered wave proportional to h-Vcpm', the superposition of these individual contributions 
forms the net scattered wave, as expected from Huygen's principle. 
Let the incident wave be a plane wave 



0,„ = e^'^■^ (A.8) 

and the surface Si be the portion of the x = plane in the interval —yo < y < 
yo- That is, it is a strip which is infinite in the z-direction. Take the wavevector 
of the incident wave to be k = k{cos9, — sm9,0) and the observation point to be 
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x' = (—6, 0, 0), that is, at a distance b from the mirror, as illustrated in Fig. |T5|. We 
can now write the scattered wave as 



0.(x') ^ ^ r ^.e---^ r (A.9) 



fx') 



27r y-j;o ^ V-oo Vy' + + 6' 

The z-integration may be performed explicitly, with the result 

'-kcosO r dye-'''''''^yH!^^\k^y^ + b^), (A.IO) 
2 "'-j/o 

where ifg^'' is the Hankel function of the first kind. In general, the remaining integral 
cannot be evaluated explicitly. However, in the case of an infinite plane mirror, 
?/o ^ C)0, it can be evaluated, with the result 

0,(x') = -e*^^'"^^^ = -e*'^'-''' . (A.ll) 

This is just the result predicted by geometric optics. In the special case of an infinite 
plane mirror, geometric optics gives the exact result. In the case of a finite mirror, 
Eq. ( A.ll| ) is still a good approximation in the high frequency limit. In this limit, one 
may evaluate the integral in Eq. ( |A.1CI| ) using the stationary phase approximation. 
So long as the classical path intersects the mirror in the interval —yo < y < yo, then 
there is one point of stationary phase within the range of integration. This intersection 
occurs if —yo < bcosO < yo. The result of the stationary phase approximation is just 
Eq. (|A.11| ), reflecting the fact that geometric optics is a good approximation in the 
high frequency limit. 

Now we wish to give a quantitative estimate of the accuracy of the approximation. 
Let 

A^, = hcos9( f r) dye-'^''^^yH^o\k\ly^ + h^)- (A.12) 

2, \J—oo Jyo / 

This is the difference between the stationary phase (geometric optics) approximation 
and the exact result for a finite mirror, so |A0s| = |A(/)s 1/10^1 is a fractional measure 
of the accuracy of the approximation. 

In the high frequency limit that k ^ b, we can use the large argument form for 

fr(l) 

//«(.) /A e^(-f), (A.13) 
V nz 

and write 

Jt^cos^c-^-/^ rrfy(y2 + 62)-3e-^'=(^-^^-v/^) + (y^ ^ -y^) , (A.14) 

V 27r Jyo 

The above integral still cannot be evaluated explicitly, but we can estimate it as being 
of order l/{ky^/y^) when k <^ y^. With this estimate, we find that 

\A<jy,\^o(-y , (A.15) 
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FIG. 15: A reflecting strip lies in the interval — yo < U < Vo ^^'^ is infinite in the z-direction, 
normal to the page. A plane wave (j^in is incident at an angle of and creates the reflected 
wave (ps at the point x', located a distance b from the strip. 

where A is the wavelength of the incident wave. 
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